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SUMMARY 

The  authm-’s  untruncated  orbital  thetxy  for  perturbations  due  to  the  Earth’s  zonal 
harmonics,  previously  developed  to  cover  all  terms  associated  with  /2^  and  // ,  is  being 
extended  to  cover  the  secular  and  long-periodic  effects  associated  with  and  If/2 

and  //  (/  >  2)  are  regarded  as  first-order  and  second-mtier  respectively,  this  means  that 
fcmnal  third-order  errors  will  no  longer  build  up  to  second  order  over  a  timescale  of  order 
up  to  I//2  in  angular  measure.  The  extended  theory  has  been  partially  checked  out  fm*  an 
Earth  model  involving  just  J2  and  Jz  . 
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AAS  91-466 


ORBIT  PERTURBATIONS  DUE  TO  AN 
AXI-SYMMETRIC  GRAVITATIONAL  FIELD 
ANALYZED  OVER  EXTENDED  PERIODS  OF  TIME 

R.  H.  Goodng* 

The  author's  untruncated  orbital  theory  for  perturbations  (kie  to  the 
Earth’s  zonal  harmonics,  previously  developed  to  cover  all  terms 
assodaledwith  and  J/.  is  bein^axtended  to  cover  the  secular  and 
long-perkKic  effects  assodsfted  with  and  J^Ji.  »J2  mi  Ji  (/>2) 
are  regarded  as  first-order  and  second-order  respective^,  this  means  that 
formal  third-order  errors  wil  no  lortger  build  up  to  second  order  over  a 
timescale  of  order  up  to  in  angular  measure.  The  extended  theory 

has  been  partialy  checked  out  for  an  Earth  model  involving  just  J2  and 
J3. 

INTRODUCTION 

The  author’s  eaiiier  wcnk*  cm  eccentiici^-untiuncated  penurbatkxis  due  to  die  zonal 
harmonics  J2  and  7$  was  subsequently^  geneandized  to  7/  and  present^  at  the  1989 
Astrodynamics  Conference,  more  detailed  verskms^*^  of  R^.  1  and  2  also  being  available. 
Refs.  3  and  4  were  conceived  as  Parts  1  and  2  of  a  trilogy.  Part  3  being  envisaged  as  the 
full  account  of  die  work  that  is  introduced  in  die  present  pqrer. 

The  main  feature  of  die  audwr’siqjfnoadi  has  been  die  particular  way  in  which  SC, 
the  short-period  perturbation  in  die  oscid^g  element  C  (^nerkfordieos^dienaents 
^e,i,  Adi  and  AOf  is  separated  frnn  the  mean  element  where  ^  C  ^ 

C,  which  have  Icmg-periodic  as  well  as  secular  variation,  are  defined  in  such  a  way  diat 
the  remnant  ^  can  be  combined  into  corxqiact  and  non-singolar  penurbatioiis  (Sr,Sb,Sw) 
in  a  particular  system  spherical  polar  coonfinat^C  with  r  s  T+Sr  eic),tQgedierwidi 
the  corresponding  patutbidions  in  velocity  (&',Sb,8wy,  analysis  for  ^  general  7/  leads^ 
to  unique  definitions  for  the  quari-constantsioqdicidy  present  in  the  (,  the ‘constant’ for 
a  being  felicitously  such  as  to  make  dw  mean  semi-ini^  axis  exacdy  constant  and  identi¬ 
cal  widi  die  quantity  a'  defined*'^  sudidiat  -iifla'  is  the  energy  int^raL  Thus  die 
thewyof  Reft.  1  aiid2naay  beregadedasha^gtwocornpotientt;  one  provides  an 
algorithm  for  the  variation  die  thereby  constituting  an  integration  c^ft^ulae 
derived  fm- die  ^  ;  the  odierconaponent  provides  ftamulae  for  A',fe,(ftv(togedigwidi 
Sr,Sb,Sw).  A  recent  piqiet^  (presented  at  die  1990  Astrodynamics  Conference)  showed 
how  the  theory  can  be  given  logical  corifAttiiess  by  three  patticular  extensions:  first 
(sonoewhat  academical^)  die  formidae  can  be  genertiHaed  to  cover  die  tiftsend  harmonic 
7bM.  so  long  as  the  Earth’s  rotation  is  ignored;  second,  die  formulae  for  ooncqitoal  7/ 
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(and  by  futiier  extension,  /bn )  widi  negative  /  can  be  obtained  (widi  an  qiplication  to 
luniscdar  perturbatitms);  third,  die  Fonran-imptonented  algcridim  can  be  extended  to  cover 
hyperbolic  arbhs. 

Unlike  Ref.  2,  die  present  paper  is  mainiyc(»cexiiedwidi  the  first  cooqpooent  of  die 
theoiy-the  ^variadon.  In  reg^  to  the  sectmd  oxiqioi^t,  it  is  noted  that  the  ji^,and 
hence  also  A*,  A  and  &v ,  are  essentially  Poisson  series  in  the  generic  angle  kcf+jv, 
where  k  and  j  are  integers,  oi' » to-  V2%  and  v  is  true  anomaly;  they  are  obtained  by 
integration  of  Lagrange’s  planetary  equations  with  V  as  integration  variable,  terms  with 
j  =  0  being  introduced  (as  tequii^)  as  the  integratitm ‘ctmstants’.  Tennswith  j  =  0 
die  integratkxi  define  die  (  ,  diou^  they  actually  arise  as  ,  diere  being  no 
fundamental  distinction  between  long-periodic  variation  (k#  0),  occurring  fix’  all  C  other 
than  a,  and  secular  variati(Mi(i(;  a  0),  occurring  for  0,(0  tmd  M  (»ly  (since  the  planetary 
equations involve_sin(ko)'+yv) fin-  a,e  and  t,but  cos(k(o' +  jv)  for  Ood  and  M). 

To  propagate  the  C>^>ccial  action  is  required  to  cqiewiA  the  obvious  singularities 
associated  with  e  s  0  and  sin  i  =  0 ,  and  the  more  esoteric  ones  associated  with 
sin^i  s  0.8  (the  critical  inclinations);  but^  this  actirm  can  effectively  be  taken  at  the 
subroutine-coding  level  and  does  not  involve  the  abandonment  of  classical  elements. 

The  existing  theory  has  been  regarded  as  ’fimnally  cooqilete  to  second  ordm’,  relative 
to  /2>  on  the  basis  that  for  />2  and  thtt  all  terms  with  coefiticient  are 

includedt>3.  The  implication  is  that  die  errors  are  of  third  order,  but  dus  is  only  true,  in 
practice,  fin- oriiital  arcs  limited  to  the  order  of  one  radian.  Over  a  timescale  with 
nt  3!  0(l//2).  in  fiict,  where  ji  is  the  mean  motion  (and  =  /t),  the  modelled  first-  and 
second-tmder  variations  in  ^  build  up  to  zero- and  first-cnderefftKts  respectively,  and  the 
uniriodelled(fininal)  third-order  variation  leads  to  second-order  error.  For  accurate  model¬ 
ling  over  quite  modest  periods  of  time,  therefixe.  there  is  a  need  to  incoiporate  the  secular 
and  long-periodic  pertubations  fhat  are  formally  of  diird  (xder,  Le.  terms  that  have  or 

J2J1  as  a  coefficient  The  third-order  perturbations  for  £i  and  m  would  involve  terms 
that  are  quadratic  in  r,inducedby  the  second-order  long-periodic  variation  in  e  and  i, 
and  effectively  leading  to  first-orto  effects  over  a  timescale  with  nt  ^  0(l//2)>  were  it  not 
fix-  the  inclusion^  of  die  apfxopiiate  ’induced  oorr^oneats’  within  the  seoxid-ocder  theory. 
It  is  instructive  to  tmdersouid  why  induced  terms  were  requited  fix:  and  m  but  not  for 
M:  they  (xiginate  fitxn  the  first-order  (secular)  formulae  for  O  and  ^  (see 
equation  (7)),  but  the  oone^iooding  formula  for  die  perturbation  in  M  (for  J2  (Xily^) 
haiqiens  to  be  null;  to  first  order,  in  Akx,  (where 

The  derivttion  of  the  secular  and  long-period  compooeats  of  die  third-order  variation, 
associated  with  coefiBdents  and  {for  /»3  in  particularX  is  die  subject  of  the 
present  paper,  the  technique  employed  bring  esseotiady  die  same  as  was  used^  for  all  die 
perturbations  associated  with  .  A  oompleie  set  of  fannilae  has  not  been  obtained, 
however,  and  success  has  been  limiied  for  reasons  that  win  emerge. 

As  in  Ref.  2,  it  is  convenient  to  conclude  die  introduction  widi  some  lenoarls  on 
notation.  The  quasi-riements  yr,  p  and  L  are  nsefid,  as  before,  wine  1^  di-f  cQ  (for 
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cs  cos  I,  and  we  also  write  faerinO,  p»  ff-f  (^iriieie  a  is  such  that  M^a+i, 

J  beiagshc  thandfor  Jndr.and  i.e^  and  Jii+  qy .  For  conciseness  we 
also  write  /  for  5^ ,  g  for  1  -  j/  and  A  for  1  -jf  .  The  lascot  p(=  acp^  and 
ii(sv  +  <o)  is  normal,  but  we  also  define  y^^pc^,  P=p/ri^  1  -t-ecosv)  and 
u'-u-jK.  As  dwp^ier  is  almost  entirely  omcemedwidi  mean  elonents,  we  again  omit 
the  bars  on  the  right-hand  sides  of  equatkms  whatever  tins  is  possible  without  causing 
cmifusimi.  Finally,  we  follow  Ref.  2  in  the  use  of  dj  and  ^  for  the  cosiiM  and  sine  of 
p  +  ku' ,  smnetimes  sumnessing  die  supeifix,  and  we  introduce  71  and  q  as  concise 
versi(»s  of  and  54  ,  that  is,  as  cos  kto'  and  sin  Am'  re^tectively.  (NB:  tiiis 
notatkm  should  not  be  confused  with  the  other  meanings  of  7  and  a.) 

FURTHER  BACKGROUND 

Following  Ref.  2  we  express  die  potential  due  to  Ji  as  Ui  -  with 

Uf  ^-Ki (pip)  Alt cos  Am'  ,  (1) 

where  Ki  =  Ji  (Rlp^ ,  R  being  the  Earth’s  equatorial  radius,  and  is  an  inclination 
function  (which  in  Ref.  2  subsumed  AT/ ,  but  sqtaration  is  preferable  now);  hoe  /  ^  2  but 
it  will  often  be  assumed  that  /^2;  further,  |■i^2  will  (as  in  previous  papers)  be 
expressed  as  K.  The  index  A  has  the  same  parity  as  /  and  0  ^  A  ^ /,  but  two  remarks 
apply  here.  First,  though  Aft  is  tmly  defined  when  /  and  A  have  the  same  parity,  a 
parallel  functitm,  A/i^ ,  was  found  useful  when  die  shot-periodic  quantity  Sb  was 
analyzed,  where  k  has  die  opposite  parity  to  /;  Ref.  2  relates  boA  the  Aft  and  the  Aix 
to  a  single  family  of  inclinatkm  functkxis,  the  A^  (0 « where  A  now  has  either  parity. 
Second,  thou^  the  restriction  to  noi-negative  A  may  seem  natural  (and  avoids  the 
duplication  of  terms  in  various  fonnulae),  it  leads  to  an  unnatural  facto  2  in  the 
definition  of  Aft  when  A  >  0 ,  and  in  the  extension  to  tesseral  harmoiics^  a  distinction 
bem  een  positive  and  negative  values  is  essential;  in  die  presem  analysis  we  effectively 
allow  for  negative  A ,  with  Aft  »  Ai^^ ,  but  it  is  only  fo  A  =  -1  that  such  quantities 
figure  in  the  final  results.  In  addition  to  Aft ,  we  n^  Aft  and  Aik  >  tbc  derivatives  with 
respectto  1 ,  though  it  is  moe  convenient  to  replace  Aft  by  Aik  =  <rhAik.  We  will 
often  suppress  the  first  suffix  in  all  these  functions,  writing  Ak  etc. 

To  facilitate  the  substitution  of  die  approfxiate  partial  derivative  of  in  each 

planetary  equation,  we  introduce  the  eccoitticity  functions  Bff ,  defined  by 

pi-i  m  (l+eco8v)^*l  ■  J!fi^cosp,  (2) 

where  overwhelming  advantage  derives  fiom  allowing  n^ative  j ,  widi  Bij  »  R/  .y .  The 
series  effectively  runs  fiom  to  -(■••,widi  By*0  when  l/l^/.  We  also  require  the 
derivatives.  By  ,  and  (as  with  the  A  functions)  we  will  (dten  suppress  die  first  suffix;  a 
useful  formula  for  not  given  in  previous  pqiers,  is 

By  -  ^(/-l)(B#.i,>fI+^M,/-l).  (3) 

In  viewof  thefiequentretpiiiementfor  R/.i,y  intfaeanalTTsisfor //  we  write  diis  in 
abbreviated  form  as  R.j . 
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The  By  are  defined  with  in  equation  (2),  dKxigh  it  is  that  occurs  in  (1), 
in  anticipation  of  the  chai^  of  vaiiaUe  in  die  planetary  equations,  the  relation  betwem  the 
variables  being 

^  =  nr^P^  +  OiJiUi).  (4) 

where  l>2  since  die  absence  a  tenn  in  J2  has  already  been  remaiked.  Afierthe 
change  oi  variable,  the  integration  for  a  fbnnal  secmd-orto  solutkm  proceeds  entirely 
straightforwanily,as£Brasthetennsin  //  in  die  sohitioo  are  cmcen^  (diese  terms  being 
first-order  in  the  second-order  quantity  Jd,  and  fiurly  strai^tforwardly  (though  laborious¬ 
ly)  for  the  terms  in  .  The  short-periodfe  terms  constitute  die  SC,  fiom  tdiich  Sr,  Sb 
Sw  are  derived,  and  are  only  of  interest  here  in  duu  th^  are  needed  to  feed  back  for 
the  third-order  solution.  The  remaining  terms  constitute  the  secular  and  k»g-periodk; 
perturba^ns,  and  in  Ref.  2  their  ccnnbinatkm  was  expressed  via  formulae  for  the  rates  of 
change  ^  .  However,  it  must  not  be  overiooked  diat  the ‘averaged*  parts  of  die  planetary 
equations,  afterjhe  change  of  variable,  cotreqiood  to  integrated  quantities  that  are  linear  in 
V  rather  than  Af,  so  that  there  are  additknial  perturtotions  terms,  slum-periodic  in  nature, 
given  by  (  C  /^)(^  -  The  resulting  difficulty  was  handled  in  Refs.  1  and  3  by  the 
concept  of  *semi-mean’  elements,  C  •  ^  obtained  from  die  C  incorporation 

of  these  additimial  terms. 

Three  ai  the  general  C  Uc  fonoudae  (associated  with  are  rqieated  frmn  Refs.  2 
and  4  (with  bars  on  the  ri|^t-hand  sides  omitted  fw  the  first  time).  These  formulae,  in 
particular,  are  required  in  die  sequel,  and  die  first  two  formulae  may  be  expressed  via  a 
single  equation.  Thus 

*  kKin<p-c  AkBkOk  (5) 

Rtk~KinrWAkB’kyk\  (6) 


for  L  ,  as  given  in  Ref.  2,  we  rqilace  er^q^Bk  , in  (6),  by  -(2/-l)qR*. 

The  cmiqileie  secular  and  long-periodic  stdution  associated  widi  (from  Refs.  1 

and  3)  is  also  repeated,  as  die  fonnulae  invtdved  are  qiecial  cases  the  formulae  to  be 
associated  with  7^/ :  they  are  ‘qiedal*,  rather  than  ‘particular’,  as  the  genoal  formulae 
give  double  the  correct  results  on  setting  /  s  2 ,  since  (as  will  become  clearer)  diey  cover 
jfJl  as  well  as  Z^// !  As  the  genenl  formulae  will  be  eqiressed  in  terms  of  KKi ,  rather 
than  K2K1 ,  we  express  the  formulae,  somewhat  artificially,  in  terms  of  KK2  (where 

K  si  Ar2);  the  weU-known  first-order  formulae  may  be  written 

^2j0  •  -Knc  and  ©2,0  =  2 iC n g .  (7) 

As  the  notation  X  21k  will  be  used  to  refer  to  die  tenns  in  or  ok  that  emerge  in  the 

solotionfior  7^/ ,  we  express  the  secular  part  of  the  known  solution  for  via  die  three 

formulae 
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f52AO  =-T6^^2nc[4(l-/)  +  e2(4  +  5/)], 

(8) 

&2X0  =  ^  «  [2(64 -18Qr+ 95/2) +  «2(56- 3^- 45/2)1 

(9) 

and 

^2.2.0  =-^Arkr2«93(8-8/-5/2). 

(10) 

It  is  possible  to  dispense  with  equaticm  (10)  by  mcxlifying  Kqiler’s  third  law  for. mean 
elements^ A  but  the  ‘additional  perturbatitm*  (term  inV-M)  resulting  from  M-y  -ip  will 
still  be  required.  The  formulae  from  which  the  long-poiod  perturbatkMis  may  be  derived 
are 

?Z^^.2,2  =  »2.2.2  =  ^KK2ne^q‘^cf(14-15f)<y2. 

(11) 

152.2,2  =  5-  KK2  n^cO-  lSf)y2  , 

(12) 

¥2.2:1  =  j^KK2n(2  +  5e2)/(i4  -  15/)y2 

(13) 

and 

^2,2,2  =  ^KK2ne^qJ{l4-l5f)Y2- 

(14) 

The  short-period  perturbations  will  not  be  ^ven  again  here,  but  two  errors  in  Ref.  1  should 
be  noted:  the  factors  r  and  J  were  omitted  from  its  equation  (43),  the  formula  fOT  Sb ; 
and  the  tenns  with  oveiaU  factor  4/  in  its  equation  (44),  the  formula  for  &v ,  should  all 
be  negated  (the  wrnig  sign  was  attached). 

To  amplify  the  remarks  in  die  Introduction  concernii^  ‘induced  terms’,  it  can  now  be 
stated  that  these  terms  arise  fiom  die  variatimi  of  ?  and  i  (via  if ,  r  and  ^)in 
equatiois  (7),  the  equations  for  this  variation  being  (S)  and  (11). 

J2J1  PERTURBATIONS 

The  analysis  for  each  element,  C  •  ^  separated  into  twp  distinct  pans  or  ‘halves’, 
which  lead  to  die  two  conqionents  of  the  final  fonnula  for  ^  2tt  •  In  the  first  half 
we  substitute  the  first-order  solutimi  associated  widi  J2  >  involving  terms  in  ^  ,  in  the 
planetary  equation  associated  with  // ,  whilst  in  the  second  half  we  substitute  die  existing 
solution  associated  widi  //  in  die  equation  associated  widi  J2  •  (The  analysis  also  covers 
/  -  2 ,  though  there  is  then  only  one  set  of  terms,  which  is  why  this  case  is  special  radier 
than  particular.)  The  index  k  is  associated  with  //  ,  not  /2 1  so  for  /2  h  is  best  not  to 
dunk  of  any  separatitm  between  ‘k^O’  and  *k  2’ ;  but  it  is  also  best  not  to  think  of  k 
as  irrevocably  associated  with  ,  in  the  second  half-analysis,  as  k  has  to  be  redefined 
during  the  develo|mient 

We  start  by  considering  genenlaqiectstrfjust  the  first  half-analysis.  Infsinciple, 
every  element  has  a  first-order  rgaesentarion,  such  as  di2,s(diere  AQ  consists 

of  ^  and  a  secular  term,  both  having  as  a  factor,  and  these  rqnesentatkms  can  be 
substituted  in  the  planetary  equation  for  (  due  to  Ji .  Because  cf  the  axial  symmetry. 
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howevor,  A  is  actually  the  one  element  that  does  not  occur  in  the  fonnula  for  ^  and  for 
winch  substitution  is  therefore  not  requited;  thus  (udy  for  n>  is  Ae  secular  compcuient 
relevant,  and  this  is  taken  account  of  sqMoately,  as  we  shall  see,  so  we  can  revert  to  the 
notatkm  typified  SQ.  The  onitribution  to  C  2J3t  arising  fixnn  SO  we  conceptually 
denote  by  Da ,  but  die  form  taken  by  the  planetary  equations  rrudces  it  preferaUe  to  work 
with  contributitms  Dp,  Dp,  Du,  Di,  De  and  Dy ,  derived  from  the  perturbations  Sp  etc 
that  are  oxnbinations  of  the  (first-order)  perturbations  in  the  elements  (so  that  'De  is  not 
die  same  as  the  />«  that  wrxild  arise  if  we  stayed  entirely  with  the  elements  themselves). 
When  ^  is  1  or  we  do  not  need  Dg  and  Dy ,  and  the  other  fimr  D  amtributions  are 


based  on  the  following  formulae  taken  direcdy  fiom  Ref.  3: 

dp  =  -jKp  \f(!eC\  +  3Cj  +  ^eC\)  ~  2h] ,  (15) 

SP  =  -j^^P  [f(3eCj  +  lOCj  + 1  leC?,)  -  4A(e  cos  v  +  3)] .  (16) 

Su  =  [2e(l  -/)Sj  +  (6  -  7/)Sg  +  2e(3  -  5/)5?,  +  4e(5  -  sin  v]  (17) 

a  =  -jKcs(eCj  +  3Ci  +  3eC^l-3).  (18) 


The  ‘constant’  in  equation  (15)  is  2h ,  rather  than  3/  (which  mi^t  be  expected  in  view  of 
(18)),  because  the  constants  in  Sa,Se  and  a  are  such  that  Sy=Yky(l-3/)  as  opposed 
to  being  identically  zero  (where  y=  pc^. 

In  regard  to  the  secular  component  of  Aa>,  this  makes  no  direct  contribution  to 
^  2/fc  •  since  its  effect  is  taken  into  account  at  the  first-order  level  by  the  evaluation  of 
short-period  perturbations  after  the  mean  elements  have  been  propagated,  so  that  a>  has  its 
value  at  time  t .  There  is  an  uuUrea  contribution  that  must  not  be  overlocdred,  however, 
arising  from  die  expression  for  the  ‘constant’  SCik(c)  •  tiie  ccmiplete  set  of  these  expressions 
being  givoi  in  Ref.  2.  Allowance  has  to  be  made  for  the  fact  that  SCik(c)  is  not  really 
constant,  specifically  (when  k  ^  0)  because  of  die  variation  of  kaf.  To  counteract  the 
effect  of  this  variation  (at  third  order)  we  incorporate  a  contribution  to  C  in  which  Yk 
or  Ok  (in  «  repiaced  by  or  -k  ®2,0  W  where  ®2.0  is 

given  by  equation  (7).  Essentially  the  same  principle  would  apply  to  the  short-period 
perturbations  in  a  complete  third-onler  analysis,  widi  die  terms  in  or  Sj  (occurring 
indie  ^  due  to  7/ )  leading  to  terms  in  k&2,oS^  or  -k&2,oC^  reflectively;  this 
principle  was  used,  with  a  detailed  explanation,  in  die  analysis  of  Ref.  3. 

In  the  second  half-analysis  the  rUes  of  J2  and  7/  ate  in  principle  reversed,  but 
expressions  fix’ the  general  Sp  (more  jxecisely  otc  are  not  so  simple  as  equations 
(1^-(18),  owing  in  part  to  the  particiilarfinn  taken  by  eadi  constant;  eadi  otpresskm  may 
be  ^v^  from  the  equations  of  Ref.  2,  however.  Tte  expression  fix*  Sp  is 
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Sp  =  -Ki  AikpiU  I,Q'  +  kyiBjCj  +  [(l+iyBk-il~l)B..k]C.k  }.  (19) 

TTic  expressicm  for  S*  may  thai  be  derived  from  its  identity  with  {Plp)(iSp  -  PSr) ,  where 
Sr{=^  is  given  by  equation  (61)  of  Ref.  2.  Similarly,  Su  is  given  by  &v-cSQ, 
where  &v  is  given  by  equadcm  (78)  Ref.  2,  togedier  widi  the  special  terms  given  by 
(8S)-(88),  and  Si2  may  be  derived  frmn  equadtms  (38)  and  (84). 

Finally  (because  there  has  been  no  time  to  take  the  general  analysis  beyond 
&  may  similarly  be  derived  from  equations  (49)  and  (83)  of  Ref.  2  the  result  being 

a  =  -Kicr^^kAik  'Lii-^ky^Bfj-^A'ikBkC.kY  (20) 

As  in  the  frrst  half-analysis,  we  must  also  ctxrect  for  the  fact  that  the  first-order  ‘constants’ 
are  not  really  constant,  because  of  effects  arising  from  the  long-period  variation  in  e,i  and 
m,  as  well  as  the  secular  variation  in  o).  The  relevant ‘constants’ may  be  taken  from 
equations  (148),  (149),  (151),  (153)  and  (164)  of  Ref.  3,  being  -^Ke(9fY2-20h)  for 
e,  jKcs  for  i,  zero  for  Q,  for  (O  ,and  ~jKqf02  for  M .  There  is  an 

important  further  source  of  terms  in  C  2tt  >  however,  which  only  appears  in  the  second 
half-analysis  and  which  we  must  deal  with  before  considering  the  analysis  for  the 
individutd  elements. 


The  terms  referred  to  stem  from  the  f>(//)  ctHtqxment  of  equaticm  (4),  which 
complicates  the  change  integration  variable  fimn  r  to  T ,  a  fuither  conqilication  being 
that  it  is  really  v  ,  the  semi-mean  equivalent  of  V ,  that  is  the  new  variable.  Now  v  is 
derived  firom  e  and  M  by  the  usual  Keplerian  procedure,  so  that 


dv 

*3"  = 


dy  de  dv  dM 


(21) 


where  dv/de  -  sin  v  (2  +  e  cos  v) ,  dvIdM  =  (1  +  e  cos  v)2  .  (22) 

—  -a  «-2  —  .  -1 

Also,  C  =  V  ^  C  .whilst  e  and  A#  (residualto  n' )  are  given  (for  the  assumed 
restrictimi  to  Uik )  by  equations  (5)  and  (6),  so  (21)  leads  to 

^  ^  nq-^p2{l-Ki  AutrHk  Bik^nvi\+P)ak-B’,k  P'^Tki)  ■  (23) 

It  is  the  terms  in  AT/  tiiat  specify  tiie  <?(//)  c<mq>on«it  of  equatimi  (4).  We  invert  (23)  to 
obtain  dt/dv  and  thence  tiie  modified  formula  fear  diangingtiie  integration  variable.  On 
expanding  the  terms  in  AT/ ,  we  get  the  fretor  that  mu^  be  applied  to  eadi  second-half 
planetary  equatimi  as  a  source  of  additional  terms  in  (m-  The  factor  may  be  expressed  as 
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-  7  *  Kie-^Ak  [(eB'k  +  kB0ieC^^2  +  4C^+l)  +  2Bi  (2  +  c2)C^ 

+  (cB;  -  kBk)  (4C.*-i  +  «C^-2)] .  (24) 

When  ^  is  X3  or  n>,  thoe  is  a  final  source  of  tenns,  (xmsideratkMi  of  which  will  be 
postpcmed  until  we  Gtxne  to  the  sectHid  half  of  die  analysis  for  i2. 

We  now  proceed  to  the  analysis  for  individual  °one  being  required  for  a .  We 
start  with  i ,  for  which  the  analysis  is  die  most  straightforward  and  fimn  which  the  result 
for  e  can  also  be  derived.  At  the  time  of  writing,  the  only  other  analysis  completed  is  for 

a. 

Analysis  for  inclination 

The  planetary  equation  for  the  first  half  of  die  analysis  is 

difdt  =  kKi  nif^csr^  Auc  ,  (25) 

where  the  right-hand  side  is  exact  as  it  stands,  the  elements  being  osculating  not  mean. 
Since  nq-^  =  (/i/p^)  >  the  right-hand  side  is  a  fiinction  of  p,  P,  u  and  i  only,  with  p 
implicit  in  AT/  as  well  nqr^ .  Hence  we  obtain  the  contributions  to  t  denoted  by 
Dp,  Dp,  Du  and  Dj  and  expressed  in  terms  of  Sp  etc  firom  equations  (1S)-(18).  In  each 
contribution  (after  differentiating  (25)  with  reflect  to  p,P,u  or  i  as  appropriate)  we 
effectively  replace  by  LBy  mi  making  use  of  equatimis  (2)  and  (4). 

Thus  Dp  is  given  by  the  terms  of  -k  (/  AT/  n p'i  c  A*  (LBjSj)  Sp  that  arc  free  of 
short-periodic  variation,  that  is,  by 

Dp  =  ^k(2l  +  3)  KKi  c  sr^  A*  {/(cBt+s  +  3Bt+2  +  3eBi:+i)oi+2 

-AhBkOk  +/(3eBifc.i +3B*.2  +  eB*.3)o]k-2}  .  (26) 


on  completion  of  the  algebra. 

Formulae  for  Dp ,  Du  and  Di  may  be  obtained  in  the  same  way  firom,  respectively, 
the  terms  of  k(l  +l)Ki  ncsr^  P"^  A*  (£Sy  5y)  SP ,  k^Ki  n  c  r^A*  (LBy  Cj)  Su  and 
kKi  n  (A'k  -/-I  Ak)(X  BjSj)Si ,  but  they  are  not  listed  here.  (The  coefficient  of  Si  takes 
this  form  because  the  derivative  of  cr^At  is  cr^A’k  -f’^A^.)  Finally,  the  effect  of 
&lk(c}  is  available  at  once  firom  equation  (20),  on  inieqireting  C.A  as  yt  and  rqdacing  it 
by  IkKngak  in  accord  with  die  rationale  that  has  bwn  noted. 

We  now  have  to  combine  five  results,  of  which  one  is  as  given  by  equation  (26). 

We  seek  a  fiuial  formula  with  0|  common  to  all  terms,  and  dim  is  no  inherent  difficulQr  in 
diis  since  the  four  contributing  formulae  effectively  tqiplyfior  off  k.  Thus  we  can  redefine 
k  (as  forecast)  such  that  die  quantity  kAkBMOk*!  in  (26),  for  exanqite,  is  refdaced  by 
(k  -  2)  Ak-TfiMOk  ■  At  this  point  we  can  conveniendy  introduce  some  notation  that  will 
simpl^  the  presentation  of  results  generally,  not  just  for  the  element  i. 
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Since  tcnm  (in  each  half-analysis,  and  far  all  elements)  arise  diat  involve  Aik+2fijt+i, 
^*+2^*.  and  Ait-2Bib.i,  we  denote 

(arbitrary)  coefficients  of  these  nine  quantities  by  X+o,  Ao+,  Aoo,  V,  X^, 
and  A.  ^respectively.  We  also  require,  for  some  elements,  similar  coefficients  for 
^it+2^Jk+i  cte  and  for  Al+^k+\  etc,  and  we  denote  these  by  /!++  etc  and  v++  etc , 
respectively.  Finally,  we  may  need  the  coefficient  of  AtB.jt ,  and  we  denote  this  by  A<.). 
We  now  introduce  tlK  generic  quantity  (2,  defined  by 

Q  -  ^-f+  ^k+lBk+l  +  —  +  /t++  ^k+2^k+l  +  —  +  A]f+2^*+l  +  —  +  Mr)  ^kB-Jc ,  (27) 

to  permit  formulae  to  be  expressed  nudnly  in  terms  of  the  particular  Q  that  is  ^qnopriate. 

In  terms  of  the  foregmng  notation  we  shall  eventually  be  able  to  write  (with  both 
halves  of  the  analysis  covered)  the  formula  we  seek  as 

i  2lk  -  /V A,  ii  c  r*  Qi  Ok .  (28) 

For  the  first  half  only,  on  combiiting  equation  (26)  with  four  other  results,  we  require 

(A++,A+o,A+-)  =  -  k-l  (ife  +  2)  {e[€{k  +  1)  -/(/  +  \Qk  +  27)], 

[6(k  +  1)  -  fill  +  7k  +  22)],  e[2{k  +  1)  -/(/  +  2k  +  7)]}  , 
(A<y+,Aoo,Ao-)  =  2{e[2(/  +  5k  +  1)  -  3/(/  +  4k  +  1)],  2[2/  -  3  -  3//] , 

<[2(/  -  5k  +  1)  -  3/(/  -  4k  +  1)] } , 
(A_^,A.o.A..-)  =  k-l  (k  -  2)  (e[2(k  -  1)  +/(/  -  2k  +  7)], 

[6(k  -  1)  +  /(2/  -  7k  +  22)],  e[6(k  - 1)  +/(/  -  10k  +  27)] ) , 

=  -  2k-l  (k  +  2)(1  -/)  (3e,  3,  e) , 

=  [0,6(6  -  7/),0],  .)  =  -  2A-I(k  -  2)(1  -/)(c,  3,  3e) 

and  all  v  =  A(.)  =  0. 

The  factor  k"* ,  appearing  in  most  of  these  expressions,  is  eliminated  when  we  combine 
with  the  second-half  results. 

The  planetary  equation  for  the  second  half  of  the  analysis  is 

di/dt  =  K  c  s  ^2u' ,  (29) 

which  leads  to  the  contributions  Dp,  Dp,  Du  and  D/  given  by  the  short-period-fiee  terms 
of  -jKnp-^  esPS^^  ,  3K  ncsS^SP ,  2Kncs  PC^Su  and 
K  n  (I  -  2f)  P  &  ;  Sp  and  & ,  here,  are  given  by  equations  (19)  and  (20).  The 
product  F  5^  (and  similarly  P  Cj )  expands  to  j  (eSj  +  25§  +  eS?|) ,  die  expansion 
being  appropriate  to  the  ^  contribution,  as  well  as  the  other  three,  since  the  factor  P  can 
be  extracted  finom  SP .  But  Sp,  SP,  Su  and  Si  involve  terms  in  dj  and  Sy  ,  the  coef¬ 
ficients  of  which  can  be  develqred  as  combinations  of  A^  Bij,  Aik  Bi-ij  and  Ajk  Bij  (for 
particular  values  of  J  lelatedto  k),  after  which  a  term  such  as  AikBySlCj  can  be 
re-expressed  as  extraction  ad tiie  short-period-fiee  pertitm  and 
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redefinition  of  k  dienleadto  j  (since  =  for  this 

term. 


The  full  second-half  analysis  for  Dp,  Dp,  Du  and  D,-  generates  considerable  algebra, 
which  will  be  omitted,  leading  to  a  result  that  (as  for  die  fint  half)  can  be  expressed  via  the 
A’s  and  fTs.  There  would  also  be  terms  in  Ak*2B~,k+l  uxl  .  wereitnot 

thatthesearecancelledby  toms  we  are  about  to  derive.  Fcdlowing  the  general  leinaiks  m 
/y/  perturbations,  we  )iave  two  additional  sources  of  terms  to  incmporate  with  the  second- 
half  ccmtributions  to  f  2^ .  First,  die  first-order  ‘constant',  viz  jK  c  s  ,  for  i  is 
respoisible  fOTcontributiws  arising  fixxn  the  long-periodic  variaocMictf  e  and  i,  as  given 
by  equation  (5)  -  the  variation  of  n)  is  irrelevant  as  (O  does  not  appear  in  this  particular 
cemstant  These  contributions  combine  to  -  KKi  ncs'^  (1-  6fiAkBk  Ok  • 

The  other  additional  source  is  the  <me  associated  with  eiqnession  (24)  that  only  arises 
in  the  second  half-analysis.  On  multiplying  this  expression  by  die  right-hand  side  of 
equation  (29),  and  then  piddng  out  the  short-period-ftee  terms,  we  get 

iKKincsAk  [(kBk  +  eB},)  Ok-i  +  (AB* - eB^)  Ok*2]  ■ 

On  redefinititxi  of  k  and  simplification  via  equatitxi  (3),  this  leads  to  ^  (/  -  1)  KKi  necs 
(Ak+2  B.^k+i  -  Ak-2  B.,k-l)<^  •  which  involves  die  cancelling  terms  t&t  have  been  alluded 
to. 


The  second-half  values  of  the  A’s  and  /r's  are  now  available.  On  combining  them 
with  the  first-half  values,  we  have  the  final  stdution  for  T  21k » expressed  via  equations 
(28)  and  (27).  The  final  formulae  for  the  A’s  and  /i's  are  as  follows: 

(A++,A+o,A+.)  =  -  {e[6(A  +  2)  -/(/  +10*  +  21)],  [6(*  +  2) -  f(2l  +  7*  +  16)], 

c(2(*  +  2)-/(/  +  2*  +  5)])  ,  (30a) 

(Ao+,Aoo,V)  =  2{e[2(/  +  5*  +  1)  -  3/(1  +  4*  +  1)],  2I2(/  -  3)  -  3/(/  -  6)]. 

e[2(/-5*+  1)- 3/(/-4*  +  1)]} ,  (30b) 

(Ah.,A.o,A-)  =  [em  -  2)  +/(/  -  2*  +  5)],  [6(*  -  2)  +  fill  -  7*  +  16)], 

«[6(*  -  2)  +/(/  -  10*  +  21)])  ,  (30c) 

Qi++4i+o4i+-)  =  -  2(1  -/)  (3e,3,e) ,  /too  =  6(6  -  7/)  , 

= -2(l-/)(e.3.3e). 

Results  for  particular  7/  may  be  doived  by  substituting  in  equations  (30).  Fot  1  =  2 
we  have  a  special  situatkm,  as  alreiuiy  noted,  su^  that  the  correct  result,  given  by 
equation  (1 1),  is  produced  by  either  half-analysis  on  its  own  and  equations  (30)  give 
dirableth^  result  But  dus  remark  posits  an  inqroitantpmnt  that  arises  for  all  evoi  1:  the 
value  of  A/.0 ,  to  be  used  as  Ak.2  when  k  =  2,  trust  tu  doubled,  to  include  the  ‘urmatural 
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facun:  of  2’  (see  the  Introduction)  that  normally  applies  only  to  with  A;  ^  0 ;  this 
allows  for  the  lestiicticHi  of  results  to  it  2  0  a^  the  consequent  absence  of  i4/,o  in  the  idle 
of  Ajfc+2  for  k  =  -2.  (This  doubling  does  not  iq)ply  to  the  use  of  A/,o  as  A*  for 
^  =  0,  but  this  is  iirelevant  for  ^=i, since  ao  =  0 ;  the  factor  of  2  present  in  A/,2  niakes 
it  unnecessary  to  include  A/,  .2.) 

Finally,  we  substitute  for  1  =  3  in  equations  (30).  with  it  =  3  and  1 ,  so  that  (27)  and 
(28)  provide  results  for  J2  ,but  first  we  note  a  point  for  odd  /  in  general.  It  is  that  for 
^  =  1  we  require  A;,  3,  A/,  1  and  A/,  -1 ;  althou^  A/,  -i  =  A/,  1 ,  we  cannot  just  combine 
the  two,  because  the  corresponding  X ’s  will  be  quite  diffoent  As  die  pair  of  formulae 
for  J2  are  new,  we  express  them  with  and  <Ti  replaced  by  cos  3o>  and  -cos  o) , 
respectively.  Then 

i  2,3,3  =  - 1*  KK2  n  cf  (20  -  21/)  cos  3m 

1 2,3,1  =  ^KK^nec  [2(8 -  15/  +  5/2)  +  5e2/  (6 - 7/)] cos  m 

These  formulae  were  originally  derived  by  a  specific  analysis  for  >  and  have  been 
checked  out  by  inclusion  in  the  authOT’s  Fortran  test  program. 

Analysis  for  Eccentricity 

The  analysis  for  e‘2ik  would  be  ccmsiderably  marc  conqilicated  than  for  f  2/;b  if  we 
proceeded  ab  initio  in  the  same  way,  but  there  is  a  much  shorter  route  that  makes  use  of  the 
exact  constancy  of  the  osculating  quantity  /(«  pc2) ,  which  is  a  measure  of  the  angular 
momentum  about  the  axis  of  symmetry.  Now 

dyfdt  =  ciq^a~2aec  e-2a  5  i  )  ,  (32) 

so  we  can  get  e2ik  from 

e2ik  = [Tm +{2pcsy^Y2ik^-  ^^3) 

It  remains,  therefore,  to  do  the  two  half-analyses  for  the  mean  quantity  7 ,  the  variation  of 
which  stems  entirely  from  the  general  non-zero  ‘constant’  Syo^c) ,  to  be  derived  for  the 
first  half,  and  the  specific  value  jlCyil  -  3f) ,  given  after  equation  (18)  and  required  for 
the  second  half. 

From  the  expressions  for  Sauqch  &tt(c)  attti  Siut(c)  io  f^cf.  2,  we  can  derive 

Syikic)  =  -Kir  lAk[il+l)Bk-(l-l}B.j^+AkBk]7k,  (34) 

from  which  the  first-half  contribution  to  yia  is  obtained  oa  replacing  yk  by  IkKng  Ck  • 
From  the  value  \k  /(I  -  3/),  similarly,  using  equatitm  (5),  we  derive  the  sec(xid-half 
contribution  such  that  in  combination  with  the  effect  of  (34)  we  get 
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ym  ^-\kKKinr[ [4(3/ - 7) - 3/  (5/ -  19)]  At Bt 

+  \2g[A'tBt-{l- \)AtB.j,)at.  (35) 


Now  equation  (35)  can  be  interpreted  via  values  of  ^oo.  /^oo  X(.) .  If  it  were  not 
for  these,  the  required  values  of  the  A ’s  and  fi 's  could  all  be  taken  as  klendcal  with  die 
values  resulting  jfrotn  the  i  analysis  since,  in  view  of  equations  (28)  and  (33),  we  naturally 
write 

=  -j^kKKin  q^QeCft.  (36) 

Thus  the  only  values,  to  be  substituted  in  (27),  that  are  different  from  those  given  by 
equation  (30)  are 


^  =  -2[8(/-2)-3/(3/-7)l, 

/too  =  12(1-/)  and  A(-)  =  24(/-l)g. 


(37) 


For  1  =  2  we  get  the  result  (doubled)  given  by  equation  (11),  whilst  for  /  =  3  we 
obtain  the  new  formulae 


«2,3.3  =  KKj  n  sf  (20  -  21/)  cos  3g) 

^2.3.1  =  ^KK^nq^s  [2(8  -  41/  +  4(^2)  _  5^2/  (6  _  7/)]  cos  (O 

As  with  the  results  for  i ,  given  by  (31),  these  fcnmulae  were  originally  derived  by  a 
specific  analysis  and  the  Fortran  test  program  has  shown  them  to  be  cmrcct. 

Analysis  for  Nodal  Right  Ascension 

For  Q  we  proceed  as  for  /,  with  a  complication  in  the  second  half  of  the  analysis. 
The  planetary  equation  for  the  first  half  is 

dO/dt  =  -Km  q-^r^  A’lt  Cj  .  (39) 

Formulae  for  Dp,  Dp,  D^  and  Di  can  then  be  obtained  fiom  the  short-period-fiiee  terms  of 
(/  +  |)  AT/  n  p-1  r  1  a;  (I  Bfj)  Sp ,  -{l+l)Kinr^  A»-l  A't  (X  BjCj)  SP, 
kKin  s-1  A’t  (X  BjSj)  Su  and  -  AT/  n  rl  (A  *  -  c  )  (X  BjCj)  Si  .  Finally,  the 
effect  of  S(2it(c) ,  given  by  equation  (84)  of  Ref.  2,  is  -*2  ATAT/  ncr^g  AtBtJk- 

We  express  the  final  fcmnula  fix  f32/ifc  as 

=  ^KKincf'^QoYt^  (40) 

with  to  be  taken  fiom  the  tqtpropriatevetsioa  of  equation  (27);  die  qrpearance  of 
c/'i ,  rather  than  5*^ ,  reflects  the  usage  of  At  ratherdian  Aj^.  Onoombining  all  die 
first-half  contributions,  and  after  the  algebra  involved  in  the  redefinition  of  A  etc,  we  find 
we  require 
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(/i++^,/i+.)  =  {e[6(*  +  1)  -/(/  +  lOit  +  21)1,  [6(*  +  1)  -/(2/  +  7jfc  +  16)], 
e[2{k  +  1)  -/(/  +  2^  +  5)]).  ()io+./iooJio-)  =  -  2{<[2(/  +  5*  +  1)  -  3/(/  +  4*  +  1)], 
2[(2/-3)-3/(/-l)].e[2(/-5*+l)-3/(/-4b+l)]).  = 

-  {e[2(*  -  1)  +/(/  -  2*  +  5)1,  [6(*  -  1)  +/(2/  -  7*  +  16)].  e[6{k  -  1)  +  /(/  -  lOit  +  21)] ) . 

(v++,v+o,v+.)  =  2/(3e.3.e).  (Vo+.Voo.Vo.)  =  (0.-iy,0),  (v_^,v^v_)  =  2f(e,3,3e), 
XoQ  -  -24*2^  and  all  other  A  =  0. 

The  planetary  equation  for  the  second  half  of  the  analysis  is 

SQ/dt  ^  -  K  n  <r^  c  (cos  2u' +  1) ,  (41) 

and  this  leads  to  Dp,  Dp,  Du  and  Di  pven  by  the  shmt-period-free  terms  of 
jKnp-^cP(Cl+l)Sp, -3Knc  (Cg+ !)£?»,  IKncPSlSu  and 
Kn  c  P  (C^  +  1)& ,  to  be  combined  via  the  redefinition  of  it  as  usual  (laborious  detail 
(xnitted).  Since  the  first-onler  constant  is  zero  for  £2 ,  tiieie  is  ix>  aridirinnal  attribution 
associated  with  equation  (5).  We  do  have  the  source  associated  with  equatioi  (24), 
however,  at  applying  this  to  (41)  and  picking  out  the  short-period-frec  terms,  we  get 

jKKi  n  e-^c  Ak  [3e(eB'k  +  kBk)Yk-2  +  4(2  +  3e^)Bkrk  +  3c(cfii  -  ikfi*)n+2]  • 

The  preceding  expression  introduces  an  apparent  singularity,  since  the  coefficient  of 
Yk  involves  a  non-cancelling  factor  outside  the  square  braclmts.  There  is  a  final 
contribution  m  the  seomd-half  li2lk  >  however,  which  did  not  arise  for  i  or  e , 

This  final  contribution  derives  fitom  the  first-order  secular  variation  given  by  ^2,0  in 
equation  (7),  wi^  a  smilar  effect  arising  in  the  second  half-analym  when  C  -  n>  and 
deriving  from  S2.0  •  This  quantity  is  a  function  of  W  (via  ^),  7  (via  D  and  n .  The 
long-periodic  variations  of  ?  and  i  were  already  accounted  for  (as  has  been  noted)  in 
the  second-order  theory,  via  die  so-called  ‘induced  components’,  errors  that  are 
0(KKi  rfl  being  thmby  avoided.  Now  die  false  assunqition  that  71  =  n'  only  involves 
CKKf^i  nt)  effects;  these  were  not  previously  relevant  but  ate  of  precisely  the  order 
covered  in  the  present  paper. 

_  The  pcnnt  is  diat  £i2,o  originates  in  the  term_-  ATr  in  diild9 ,  and  tltis  integrates  to 

-  AT  7,  which  is  handled  as  the  secular  term  -ATr/i'r.with asa 

separate  short-periodic  ^ect  But  has  an  0(Ri)  compoaent  denying  fctmi 

equations  (5)  and  (6)  via  (21)  and  (22).  Equatitms  (5)  and  (21)  lead  to  no  long-term  effect, 
but  (6)  and  (22)  lead  to  our  requiied  ‘final  amoibution’  of  -jKKi  ne-^c  (2  +  e^) 

AkB'kYk  on  taking  the  mean  value  of  p  be  j4r3(2  «2).  In  cmnbination  with  the 
preceding  contribution,  this  effectively  changes  the  co^cient  of  B'k  Yk  (within  the  square 
brackets)  from  4(2  +  3e2)to  Se^ ,  thus  eliminating  the  af^iarent  singularity. 

It  remains  to  ooodnne  all  die  contributioos  in  die  second  half-analysis,  to  redefine  k, 
and  to  oondiine  widi  die  results  of  the  first  half-analysis.  In  tenns  of  equations  (40)  and 
(27)  we  finally  derive: 
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(^++.A+o,V)  =  2/{«(/  -%k-  21).  -  (2/  +  3*  +  16).  -  tf(3/  +  5)) .  (42a) 

(A0+.A00.A0-)  =  -  6{ef(l  -5k-  2).  [4*2  -/(5*2  +  4/ -  8)3.  ef(l  +  5k-2)),  (42b) 

(A-+.A.O.A..)  =  -  2f[e(3l  +  5).  (2/ -  3*  +  16).  ^(/  +  8* -  21)} ,  (42c) 

(/z++./i+o.M+-)  *  U[6(Jk  +  3)  -/(i  +  10*  +  33)1.  16(*  +  3)  -/(2i  +  7*  +  28)]. 

e[2(*  +  3)  -/(/  +  2*  +  9)]}  ,  (42d) 

(/io+,/ioo./A>-)  =  -  2{e[2(/  +  5*  +  1) -  3/(/  +  4*  +  1)1.  2(2/ -  3  -  3//). 

«t2(/  -  5*  +  1)  -  3f(l  -  4*  +  1)1 } ,  (42c) 

(H.+.^K„H„)  =  -  {c[2(*  -  3)  ■»■/(/  -  2*  +  9)1.  I6(*  -  3)  +/(2/  -  7*  +  28)1, 

c[6(*  -  3)  +/(/  -  10*  +  33)1 ) ,  (42f) 

{v++,v+o,v+..Vo+,Voo,Vo.,  v_^.  v^,  v„)  =*  y(3c.3.c.0.-6.0.c,3.3c)  (42g) 

and  A(.)  =  12(1  - 1)/(4  +  3c2) .  (42h) 

For  1  =  2,  equations  (42).  substituted  in  (27)  and  (40).  give  the  tesults  (doubled) 
known  from  equations  (8)  and  (12),  with  *  =  0  and  k  =  2  respectively;  for  *  =  2 ,  as 
already  noted,  we  require  the  values  of  i42,o,  and  A2jo  to  be  doubled,  whilst  for 
*  =  0  we  use  these  values  undoubled  and  do  not  require  A2,.2  etc.  For  /  =  3  ,  with  *  =  3 
and  *  -  1 ,  we  get  the  new  results  given  by 


^2.3.3 

^2.3.1 


-“I  KK^  ne^cs  (4  +  If)  sin  3o) 

^KKznecs-^  [2(8  +  9Sy -  185/2)  +  5e2/(i8 _ 35/))  sin  o) 


(43) 


Equation  (43)  conftxms  with  the  results  oi  die  qiecilic  analysis  originally  carried  out, 
but  is  unfortunately  not  validated  by  the  test  program.  There  are  effectively  seven 

numerical  coefficients  in  the  equations  aiKlconqiuter  runs  different  values  of  e  and  i 

indicate  that  five  of  them  are  correo.  The  dubious  integers  are  99  and  -185,  with  some 
evidence  that  the  correct  values  might  be  51  and  —125,  which  would  imply  an  enor,  in  the 
square  brackets  of  ^523,1  >  prt^xwtional  to  /f .  At  the  time  of  writing,  it  is  not  known 
whether  the  error  is  a  one  or  an  artefact  of  the  algorithm  used  in  the  test  program.  (The 

two  leading  coefficients  in  equations  (43),  contributing  4  sin  3©  and  16  sin  ©  ^part  from 
the  overall  factors  at  die  beginning  of  the  expressions,  could  actually  have  been  written 
down  without  any  analysis,  since  diey  are  mandated  by  the  cone^ionding  coefficients  in 
equations  (31),  to  avoid  singularity.) 

Analysis  for  Psrigss  Argumant  and  Moan  Anomaly 

The  general  analysis  for  die  last  two  eloiimts  has  not  been  completed  at  die  time  of 
writing.  Specific  results  for  have  been  derived,  but  (as  for  )  they  are  not  validated 
by  the  test  program  and  the  formulae  will  not  be  quoted.  (The  formulae  for  ^233  and 
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^23,1  >  which  j&233  snd  @23,1  follow,  involve  four  and  nine  numerical 
coefiBcients,  re^ecdvdy;  two  of  the  four,  and  three  of  the  nine,  are  certainly  correct  from 
their  mandatory  aaaodatkn  widi  the  coeffidents  i^tpearing  in  equation  (38).)  In  die  present 
ptqier  we  just  remark  on  some  aspects  of  the  general  analysis. 

The  analysis  for  die  quasi-element  yr  is  distincdy  sunpln  than  for  on  itself,  rince 
the  planetary  equatkm  is  sinqiler  for  (/duuifcn'  d>(=  cl2)  and  the  analysis  fen*  Q 
has  already  been  d(Mie.  We  must  allow  for  the  vaiiadon  of  c  (=  cos  1),  but  ^  merely 
involves  the  incoqxnation  of  conttibudons  given  by  sQ  Si ,  where  £2  as  given  by  (39) 
and  Si  by  (18)  are  appre^niate  for  the  first  half-analysis,  and  similarly  for  the  second.  The 
analysis  is  necessarily  more  laborious  than  for  1  and  £2 ,  however,  as  shown  by  die 
''lanetary  equation  for  the  first  half-analysis,  whkh  (cf  equations  (25)  and  (39))  is 

d\iridt  »  -jKinq-^AucP^*^  { (/+ +  k) +  2(/+ +  2it)e-Vcf  +2l*C^ 

+  2(1+  -  2k)e-l  dti  (/+  -  k)  ^2  )  .  (44) 


where  /+  =  /+!  (cf  equation  (39)  of  Ref.  2). 

Equation  (44)  leads  to  the  usual  quantities  denoted  by  Dp,  Dp,  Du  and  D,- ,  but  we 
now  also  need  Dg  and  Dy ;  the  former  is  associated  with  the  two  appearances  of  in 
(44),  and  the  latter  with  the  implidt  appearance  of  v  in  Cy  ,  which  becomes  explkit  (if 
j  *  0)  when  we  write  dj  in  fiilL  The  expressions  for  &  and  Sv  (see  equations  (148) 
and  (182)  of  Ref.  3)  are  matt  oxiqilicated  duui  the  fmir  quantities  givoi  by  equations 
(15)-(  18),  so  we  cannot  expect  thi^  after  redefinitkm  of  k  in  the  analysis,  all  the  coef¬ 
fidents  of  "Yk  in  can  be  related  to  just  and  Rjt.i ;  the  number  of  B’s 

required  in  each  half-analysis  in  fact  rises  from  diree  to  eleven,  but  on  combination  of  the 
two  halves  it  reduces  to  three  again.  This  avdds  die  occurrence  df  ixm-zero  p21k  ^th 
k>l . 

Analysis  for  the  sixdi  element  in  practice  invdves  the  quasi-element  L  rather  than 
M ,  and  this  in  turn  involves  die  planetary  equation  far  the  quasi-element  p  rather  than 
a ,  where  p- QY-  The  analysis  fw  p  proceeds  in  the  same  way  as  for  yr, 
allowance  being  made  for  the  variation  of  q  in  qY,  and  then  a  final  analysis  is  required 
for  the  integral  of  n ,  the  quantity  defined  such  diat  L  »  p  -t-J.  The  results  for  02lk 
and  l2ik  ,  individually,  do  not  reduce  to  zero  fca*  k^ /  + 2,  even  when  the  pairsof  half¬ 
analysis  are  oranbined,  but  on  combining  these  cranbinatitHis  we  do  get  results  for  L2ik 
in  our  standard  form. 

Now  f  *  n't  iSn  dt  and  we  may  conveniently  write  Sn  «  n'(K  fi  +Ki(^  + 
kTUT/ ^2/)  to  die  accuracy  that  concerns  us  (with  additional  terms  in  and  when  we 
comefo  ^2^  perturbations).  Exfuessing  Sa  siixBlarly,  we  require  that  , 

from  which  we  obtain 

« .-|a,  ft,  --li,,  Hu  -  (45) 


TMSp386 


18 


The  second  tarm  in  die  oqnession  for  $2/  ^  regarded  as  an  easily-evaluated  ccnec- 

don  term,  so  we  consider  only  the  main  term.  Tlius  we  require  to  evaluate,  as  the  princ^ial 

part  of  Jm  >  the  diiid-rndercoaqxment  of  -  (3n/2a)  Sa .  Changing  the  integration  vaii- 
abte  as  usual,  we  have 


dJ/dP  =  -  (3/2fl)  (^F-'^Sa,  (46) 

so  we  need  to  etqmess  the  required  compcments  of  Ba  with  as  a  factor. 

We  start  from  the  exact  equatkm 

Sa  2a’ Ki<r^  Au, ,  (47) 

all  quantities  (other  than  o')  on  the  right-hand  ade  being  osculating,  not  mean.  Whoi 
these  quantities  are  expressed  for  die  first-order  scdution,  we  can  derive  the  usual  first-half 

part  of  l2tt  <  via  the  ctHiqionents  Dp,  Dp,  Du,  Di  and  Dg ;  the  need  for  Dg  stems  from 
the  factor  q-^  in  (47).  There  is  no  difficulty  in  either  half  of  the  analysis,  but  the  details  are 
omitted  as  the  conqilementaiy  analysis  for  p2lk  has  not  been  cooqileted. 

PERTURBATIONS 

The  analysis  for  pertuibatimis  follows  the  same  principles  as  fcv  72^/perturba- 
tions.  die  sect^-ocder  scdution  for  (  being  fed  back  into  the  planetaiy  equation.  The 
analysis  is  more  tedious,  however,  and  there  is  no  cqipartunity  for  cross-check  between 
general  and  particular  results.  Atthetiitieofwriting,fotmulaehavebeenobtainedfor 
^  =  i,  e  and  h ,  but  none  has  been  successfully  validated  by  the  test  program. 

Some  revision  of  notation  is  called  for,  mch  that  the  suffix  3 ,  in  particular,  implies  a 
power  of  J2  (or  in  practice  ^  radier  than  a  value  of  / .  In  general,  we  return  to  die 
notation  of  Ref.  3,  writing  tire  known  short-periodic  cmnponent  of  the  second-mder 
solutimi  as  (fij^  = )  KCi  +^2^2,  where  Cl  C2  are  Poisstm  series  in  cos  jv,Cj 
and  Fj  (or  die  conesponding  sines),  where  (widi  bars  supinessed)  Cj  ~  cos(/v  -t-  2qi)  ^ 
7}  =  cos(/v  + 4(0);  in  the  earlier  notation  of  d^  paper,  Cj  and  /}  would  be  written  -Cy.2 
and  Cjji  .  (As  in  Ref.  3,  the  notations  Cl  C2  are  extended  to  any  quantity  that  is  a 
function  of  the  orbital  elements.)  We  seek  formulae  for  the  third-order  mean  rates  of 
change,  which  will  be  exjnessed  as  C  30  secular  component  and  C32f<vtfre 

long-periodic  component;  die  latter  notation  is  used  because  the  long-period  rate  of  change 
is  prc^iottional  to  cos  2 lb  or  sin  22b,  terms  in  4m  being  absent 

As  with  the  y^/analysis,  it  is  helpful  to  eiqness  the  planetary  equatimis  in  terms  of 
the  quantities  p,P,u,i,v  and  e ,  rather  than  die  original  ekmenis;  since  the  analysis  has 
not  been  carried  to  C*<o  C*^  p,P,u  and  i  here.  The  notation  Dp 

etc  is  no  longer  hel^pfril,  however,  as  we  effectively  require  14  basic  contribntioiH  to  C  30 
and  C  32 :  foof  associated  with  P2,  Pi,  ui  and  12 ;  four  widi  Pi^,  P^,  and  ; 
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and  six  widi  die  cross-products  p\P\  etc.  There  are  also  contributioiis  associated  with  the 
‘constants’  in  Cl  C2  >  and  with  the  OiJ'^  part  of  equation  (4). 

Because  p  and  P  occm' in  the  same  way  in  four  planetary  equatic»s(C=u  and 
C-Af  are  the  excepdtms),  namely,  via  JIT  n  ,  it  is  ctmvment  u>  set  up  a  grouping 

of  the  five  basic  contribudtms  associated  with  p2,  P2,  P^,  P^  and  p\P\ .  It  is  also 
convenient  to  introduce,  fcr  p  and  P  cxily,  a  ‘normalized  notation’  si^  diat  ^  =  p\lp 
etc;  this  is  mm-trivial  for  P\  and  P2«since  P  is  itselfa  Poisson  series,  give^  by  ^ 
equaticm  (2)  with  /  -  2 ,  by  which  Pi  and  P2  must  be  exacdy  divisible  for  Pi  and  P2 
to  have  us^iil  meanings.  The  formula  for  ^1  is  given  at  once  by  equation  (187)  of  Ref.  3, 
which  also  includes  all  die  odier  Cl  and  C2  that  we  need,  an  excqition  being  P2 .  which 
is  given  by  evaluation  of 

P2  =  «2  cos  v-ev2sinv-jvi(2eisinv  +  evicosv).  (48) 

(It  needs  38  separate  cosine  terms  to  express  each  term  of  (48),  but  most  of  these  cancel  out 
on  combination;  the  remaining  18  dien  reduce  to  13  when  Psl-t-«cosv  is  divided  out- 
cf  (51b)  below  -  and  this  cancellation  is  re^ionsible  for  the  absence  of  4  m  terms  in 
C  32  •)  We  use  the  expansions  (ctxrect  to  second  rader) 

(p/p)-’^  *  (02- 

and 

(P/pf  =  l+3KPi  +  3^  (P2  +  Pi^)  , 
from  which  it  follows  that  die  grouping  of  taros  we  require  is  given  by 

3(P2  +  ^1^)  S0i  Pi) ; 

this  leads  to  an  expression  of  the  form 

2^  {/2(173e2  re  + ...  +  327<2  rn  -  4/  (terms  in  C4 . Co) 

-t- 16 (terms in  cos2v,cosv  and  1)} .  (51b) 

There  are  nine  remaining  ‘basic  conoibutions’  to  be  covered,  for  each  C  sqiarately: 
two  lead  to  a  particular  combination  of  U2  and  ui^  ;  two  lead,  similarly,  to  a  combination 
of  12  and  ;  one  comes  from  ui  ii ;  and  the  last  four  lead  to  the  product  of  an  qipro- 
priate  combinatioa  of  wi  and  ii  with  the  fixed  combination  of  pi  and  Pi  given  by 

jPl  -3Pi  -  ^(r(5eC3  +  12C2  +  9rCi)-4A(3ecosv  +  2)}.  (52) 

In  regard  to  die  contributions  10  ^  30  and^  32  due  to  the  non-constancy  of  the 
‘constants’ in  Cl  ^  (2  •  we  start  with  die  {^1  constants,  wfaidi  were  listed  after 
equation  (20).  Each  constant  has  a  varhdion  due  lo  tfie  variation  in  W  (wliicb  itffocts  die 
siqipressed^  in  JfC)  ^  7,as^iecifiedbyeqaalk»i(ll),andcanc^ation(rfdiis 
variation  is  one  source  of  the  lequired  contribution.  Since  the  constants  in  tui  and  Afi 
(dioughnotin  fior  i3i)iiiv(dve  S ,  we  also  have  10  cover  die  effects  of  S  variation. 


(49) 

(50) 

(51a) 
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with  a  coiiq>lication  assodafed  with  die  first-ocder  coa^onmt  ttf  this  variadoii,  but  diese 
effects  are  not  required  in  the  incomfdete  analysis  df  the  present  pqier.  In  dealing  widi  the 
‘ccHistant’ part  (tf  {^2 1  on  the  other  hand,  it  is  only  the  variation  of  o  that  is  relevant  -  it 
is  dealt  with  in  exactiy  die  same  way  as  in  SQtfe)  in  die  second  half  ctf  die  analy^  the 
J2fi  perturbadcMis. 

To  get  the  contribudon  to  ^  30  and  ^  32  due  to  die  part  of  equation  (4),  we 
apply  equations  (21)  and  (22)  again.  This  time  we  use  expressions  for  e  and  M  given 
by  equations  (10),  (11),  (13)  and  (14).  From  M17  0  we  get  the  factor,  analogous  to  (24), 
given  by 

^  A:2  (8  _  8/  _  5/2)  [g2  cos  2v  +  4e  cos  V  +  (2  +  e^)] ,  (53) 

whilst  from  ?  7  and  A/2,2,2  we  get 

(14  ~  15/)[e2C2  +  4eCi  +  (2  +  e2)Col .  (54) 

When  is  /2  or  ni,  thoe  is  a  final  source  of  toms  in  ^  30  and  ^  32 ,  associated 
(just  as  in  the  second  half  of  the  ^2^/  analysis)  widi  the  interpretation  of  n  in  equations  (7). 

Analysts  for  Inclination 

Lagrange’s  planetary  equation  has  been  given  by  (29),  which  can  also  be  written  as 

di/dt  *  -  ie(ff  P^)  (pipf^  {Pipf  Pc  sSi.  (55) 

Then  the  change  of  variable,  specified  by  equation  (4),  leads  to 

dild9  =  -  Kiplp)'"^^  {Pipf  Pcs  S2.  (56) 

We  have  developed  (p/pY^^  and  {P/P via  equations  (49)-(51),  so  for  the  basic 
contributions  to  T  32  (there  is  no  T  30)  it  remains  to  develop 

cs  =  rr+^(l-2  7)11  +  [(1  -  2  f)i2  -  r  (57) 

52  =  S2  +  2KC2Ui+2k^(C2U2-hui^)-  (58) 

We  now  worir  with  equation  (56)  to  derive  all  the  cmitributions  to  i  32 ,  expressing 
each  contribution  with  the  overall  factor  (V2304)  ^cs  sin  2<o  suppressed. 

From  the  grouped  terms  rqxesented  by  equations  (51),  widi  the  factor  ~K}P  csS2 
iqiplied  and  the  shott-period-fiee  component  dim  picked  out,  we  obtain  -3(448-153^+ 
Vn9f2-).  From  the  product  of  (52)  with  IfiP  [(1  -  y)  ft  *l  +  2c  s  C2  «i]  we  obtain 
4(728  -  216y  +  \6nf2).  From  the  product  of  -4^P  with  the  combination  of 
2cs(C2U2-S2Ui\  521(1 -y)f2- 2c  and  2(1  -  V)  C2  ui  ii) ,  we  obtain 
(5408- 11724^ ••■560^^.  And  the  coodiinatioa  of  there  dote  lesulis,  representing  the 
sum  of  all  die  basic  eontributkms,  is  8(872  -  1974f 105S^). 
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The  ‘ccmstant’  in  I'l ,  viz  jc  s ,  leads  (via  t'i  op  and  i  ?  9  ?)  to  -48(14  -99f+ 
9(y2),  and  the  Co  ccraqwnent  of  die  constant  in  12 ,  viz^  c  s  (9  +  f)Co ,  leads 

(via  ©2,0)  *0  -32(36  -  41/- 5/2).  And  the  combination  01  these  two  le^ts  is 
-16(114-375/+ 2^). 


Finally,  the  product  of -KP  csS2  with  the  factor  given  by  (53)  leads  to 
-36(8  -  8/-  5/2),  there  being  no  ctmtiibutitxi  from  (54). 

On  ccanbining  all  ccmtiibutions  we  get 

732  =  5^Ar3  ne^cs  (1216  -  2361/+  1 123/2)  sin  2© .  (59) 

Analysis  for  Eccentricity 

As  in  the  analysis  for  J‘^i  perturbations,  we  derive  ?32  from  7  32  via  the  exaa 
ccNistancy  of  y(  =  pc2).  To  apply  the  equation  analogous  to  (33)  we  require  the  formula 
for  ^2  ’  ^  noted  in  Ref.  3,  y  has  no  first*  or  second-order  variation. 


The  variation  of  y  stems  entirely  from  the  non-constancy  of  the  ‘constants’  in  yi 
and  72  •  The  cwistant  in  71  is  |•7  (1  -  3/),  which  leads  (via  ?22,2  and  i 2.22)  ^ 
connibutitm 

n  e^yfiU  -  15/X5  -  9/)  sin  2q) 

to  •  The  Co  component  of  the  constant  in  72  is  -^e^yf  (63  -  8y)Co ,  which 
leads  (via  ©2.0)  to  the  contribution 

-  <27/(63  -  8y)(4  -  5/)  sin  2©  . 


So 


%2  =  -  12  -  241/  +  14(y2)  sin  2(0  .  (60) 


From  equations  (33),  (59)  and  (60)  it  now  frdlows  that 

^32  *  -  « <72/(320  _  43y  +  3/2)  sin  2©  .  (61) 


Analysis  for  Nodal  Right  Ascension 

We  are  now  kxddng  for  secular  variation,  reinesented  by  Hjo  >  as  well  as  the  long- 
period  variatkni  rqnesented  by  ^32;  The  analysis  is  very  similar  to  that  for  1,  except  for 
a  final  contribution  that  stems  frtmi  D2fi  and  the  non-identity  of  n'  and  H ,  just  as  with 
JjJl  perturbations. 

The  planetary  equation  has  been  given  by  (41)  and  can  also  be  written  as 

dO/dt  =  K(n^''^P^)(pfpy'^^{PlP)^  PciC2-l),  (62) 
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whence 

A 

dQ/d9  =  K(p/py'''^(PlF)  PciC2~l).  (63) 

For  Ae  factors  involving  pip  and  PIP  vic  can  use  the  grouped  terms  represented  by 
equadcms  (SI),  whilst  for  c  (C2  -  1)  we  require,  in  analogy  with  (57)  and  (58), 

c  =  rii-ie^(r/2+ (64) 

and 

C2-I  =  C2-1-2^  J2«l-2i?^(j2M2+ C2«1^)  .  (65) 

We  now  express  the  contributions  to  /330  with  the  overall  factor 
suppressed,  and  to  iiyi  with  the  factor  (V2304)1^R  cos2ci)  suppressed. 

From  the  grouped  terms,  with  die  factm*  P  c(C2  - 1)  applied,  we  obtain 

-[16(40  -  144/+  113/2)  +  «2(576  -  2038/-I- 1831/2)]  and  3(448  -  2816/+  2785/2)  . 

From  the  product  of  (52)  with  ^  P  [s(C2  -  l)/i  +  2c  52  «i]  we  obtain 

8(4(6  -  13/+  6f2)  +  <2(58  -  159f  +  8^]  and  -8(400  -  1107/+  711/2) 

From  the  product  of  I^P  with  the  combinaticHis  of  -2c(S2U2  +  C2Uih, 

-(C2  -  l)(si2  +jcii^)  and  2sS2  u\  ii ,  we  obtain 

-[16(24  -  41/+  5/2)  +  c2(176  _  294^+  2955/2)]  and  -(6752  -  1341^+  5907/2) . 

Hence  the  sum  of  aU  basic  contribution  yields  (with  the  same  overall  factors  assumed) 

-2(8(52  -  159/+  106/2)  +  3^2(48  -  61^+  683/2)]  and  -8(1076 -  1728/+  405/2) . 

The ‘constant’ in  is  zero,  so  there  is  no  contribution  to  1^30  or  ^32  via  C2;2.2 
or  i  2A2  •  But  the  constant  in  is  - e2  (ig  -I9f)  So ,  and  this  leads  (via  0)2,0) 
to  a  contribution  to  1^32  of  16(72  -  16^  +  95/2)  (multiplied  by  the  usual  factor). 

We  treat,  finally,  the  contributions  to  iiso  and  ^32  arising  frinn  the  product  of 
K  P  c(C2  -  1)  with  the  factors  given  by  (53)  and  (54),  together  with  the  contributions 
associated  with  the  fact  that  the  factor  n  in  equation  (7)  ought  to  be  inteqineted  as 
«'  +  3?2,2,0  +  Mf,7,7 ,  but  is  in  practice  set  to  just  n’ .  The  two  pairs  of  contributions  can 
conveniently  be  combined  since  (it  can  be  shown  that)  the  effect  of  the  second  pair  is  to 
delete  the  factor  2  -t-  e2  in  (53)  and  (54).  As  a  result  we  obtain 

-18eV(14-15/)  and  72(8- y  + 5/2). 

On  ctmobining  all  contributions,  we  derive 

^30  -  -  nr  [4(52- 155/ +  10^ +  3e2(24  -  28y  +  31Sir2)]  (66) 
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and 

^32  =  (430  -  66^  +  85/2)  cos  2®.  (67) 

There  is  a  conflict,  associated  with  singularity,  between  the  leading  coefficients  in  equations 
(59)  and  (67);  this  has  not  yet  been  resolved. 

Analysis  for  Perigee  Argument  and  Mean  Anomaly 

Hie  analysis  for  the  lat  two  elements  has  not  been  embariced  upm  at  the  time 
writing. 

DISCUSSION 

An  orbital  theon^  is  of  no  practical  value  until  it  has  been  inqilemoited  in  an  accurate 
and  efficient  conqiuter  program  for  ephemeris  generadcm.  The  audior’s  untruncated 
second-order  theoyi'^  for  an  orbit  in  an  axi-symmetric  gravitational  field,  or  any  non¬ 
rotating  Held,  had  been  fully  validated  by  such  a  computer  program,  the  main  features  of 
which  are  as  follows.  First,  a  mean-orbit-based  coordinate  system  was  utilized  as  being 
ideal  for  the  efficient  representatitMi  of  the  short-period  perturbations  free  of  singularity. 
Secondly,  singuLuity  problems  in  the  prqiagation  of  die  mean  elements  (  0  were  avoided 
by  two  expedients:  two  of  the  C  formula  are  for  C=¥  ^  ^  •  >^cr  than  eo  and  M , 
with  jD  and  stored  (rather  than  D  and  p) ;  and  die  non-singular  mean  elements 
7  sin  ,  -7  cos  and  7  are  introduced  locally  in  the  propagation  of  7  and  £2  ,  to 
ccxiqilete  the  avoidance  of  any  problem  when  7  »  0  (with  a  similar  procedure  that  pays  off 
when  ?»0).  Thirdly,  the  concqit  of  a  soni-mean  element  was  employed  as  a  device 
associated  widi  die  muisformation  of  the  integration  variable  from  r  to  T  whilst  the 
prqiagation  of  the  C  ^  required  in  terms  of  t .  Fourthly,  the  terms  in  ^  and  d> 
that  are  induced  by  the  second-Wer  terms  in  ?  and  7 ,  and  are  formally  only  of  third 
order  but  responsible  for  quadratic  variatkm  with  t ,  were  included  in  the  program. 

Hfthly,  and  to  suppl^ent  the  fourth  feature,  an  option  was  introduced  to  ‘rectify’  the 
propagation  of  the  C  h  (epoch)  to  r  by  the  use  of  intermediate  epochs  as  way- 
stations;  use  of  the  c^tirni  degrades  die  status  (k  the  diemy  from  frilly  analytic  to  semi- 
analytical,  but  die  efficacy  of  the  ‘induced’  terms  is  sudi  that  the  interval  r-  <b  has  to  be 
equivalent  to  several  hundred  orbital  revdutkms  before  there  is  any  gain  fix>m  swishing  to 
the  semi-analytical  mode.  Hnally,  the  inverse  of  die  al^mthm  that  converts  the  C 
position  and  velocify  was  made  Ughly  efficient,  and  of  essentially  unlimited  accuracy,  by 
the  use  of  a  general  iterative  inversion  procedure^  (it  fypically  gives  8-<lecimal  accuracy 
after  a  single  iteiatkm!);  this  procedure  is  at  die  heart  of  the  program’s  operation,  since  it 
permits  the  conversion  of  successive  state  vectors  (ctm^onaits  of  positiem  and  velocity),  as 
generated  by  an  independent  pure-numerical  integration,  into  successive  sets  of  mean 
elements  which  can  be  compared  with  sets  propagated  by  the  theory-based  algorithm. 

The  last  feature  above  may  be  anq^edl^  two  remaiks.  I^rst,  die  inversion 
ixocedure  is  indqiendent  of  the  algorithm  it  inverts  (apart  fimn  the  frax  that  it  calls  the 
algorithm  direedy  during  each  iteration),  so  that  modifications  to  die  algorithm  (sudi  as  die 
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adding  of  higher-order  terms)  do  not  affect  the  inversion  procedure  itself.  Secondly,  the 
existence  of  the  procedure  provides  an  intrinsic  and  definitive  answer  to  die  question: 

“Yes,  but  what  actually  are  your  mean  elements?”. 

In  view  of  the  successful  validation  of  the  second-csder  theory,  originally  fw  just  J2 
and  73  but  later^  with  the  perturbations  due  to  incmporated,  no  fundamental  difficulty 
was  anticipated  in  extending  the  theory  to  cover  the  secular  and  long-periodic  effects  of 
third  order.  (That  the  algebra  would  be  long  and  tedious  Hm  anticipated,  and  the  use  of  a 
computer-algebra  package  was  consitter^  but  appeared  to  be  impracticable.)  After 
validated  formulae  had  been  found  for  i  and  e  associated  with  the  product  , 
therefore,  it  was  assumed  that  completion  of  the  results  for  ,  the  analysis  for  ,  and 
the  generalization  from  to  would  all  be  straightforw^.  The  intention  was  that 
the  updated  computer  program  would  then  be  used  for  a  full  examinati(Mi  of  the  long-term 
accuracy  of  the  extended  theory.  Since  discrepancies,  as  yet  unexplained,  have  been  met  in 
the  analysis  for  both  JjJz  and  ,  however,  die  work  is  incranplete  and  the  present  paper 
must  be  regarded  as  no  more  than  an  interim  report 
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